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Abstract 

Structure constants of the su{N) {N odd) Lie algebras converge 
when N goes to infinity to the structure constants of the Lie alge¬ 
bra sdiff{T'^) of the group of area-preserving diffeomorphisms of a 2D 
torus. Thus Zeitlin and others hypothesized that solutions of the Eu¬ 
ler equations associated with su{N) algebras converge to solutions of 
the Euler equations of incompressible fluid dynamics on a 2D torus. 
In the paper we prove the hypothesis. Our numerical experiments 
show the Galerkin method applied to Euler equation of hydrodynam¬ 
ics is computationally more efficient in the range of time in which it 
is stable than that based on the SU(N) approximation. However, the 
latter is stable for much longer time. These numerical results agree 
with theoretical expectations. 
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02.70.Dh — Finite element and Galerkin methods 
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02.20.Tw — Inhnite dimensional Lie groups. 


1 Introduction 


In the beginning of 90’s a hnite-mode approximation of 2D ideal hy¬ 
drodynamics was introduced (|]l|, P]). The approximation is based on 

convergence of the structure constants of the Lie algebras su{N) (iV-odd), 
when iV —>• cx) to the structure constants of the Lie algebra sdiff{T‘^). The 
Lie algebra sdijf{T‘^) consists of divergenceless nonconstant vector helds on 
a 2D torus. The vector helds are inhnitesimal generators of the group of 
area-preserving diffeomorphisms of the 2D torus. Incidentally, the group can 
be also identihed as the group of symplectic diffeomorphisms, for the natural 
symplectic structure of the torus. The group is the basic symmetry group for 
incompressible ideal hydrodynamics on a 2D torus. In the physical model 
the huid is maintained in a region of a plane, appropriate boundary condi¬ 
tions being satished. The approximation of sdijfiT"^) by su{N) is used quite 
widely both in hydrodynamics and in string theory § g, 1,1,1,1,0, |]. 
There are also other hnite-mode approximations of the Euler and similar 
equations related to symplectic geometry i, in, 0,0. 

In m the approach to approximate sdiff{T‘^) is formalized and a weak 
uniqueness theorem is proved. Unfortunately, the weak uniqueness does not 
lead to an isomorphism of the limit algebras. Indeed, there are examples 
of various limit algebras which have their structure constants equal to the 
structure constants of some inhnite dimensional algebras and despite this 
equality the limit algebras are not isomorphic to the inhnite dimensional 
ones. As shown in [0, the ^im su{N) is neither isomorphic to sdiff{T‘^) nor 
to sdif f{S^), where sdiff^S"^) consists of divergenceless vector helds on the 
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2D sphere 3“^. This approximation is still not well understood and additional 
work should be done towards its satisfactory formulation. 

The plan of the paper is as follows. 

In the second section we discuss the issues connected with approximat¬ 
ing the Lie algebra sdiffiT"^) by the Lie algebras su{N). We point out 
weaknesses of this approximation. 

In the third section we remind how the Euler equations are being asso¬ 
ciated with arbitrary Lie algebras. Then we specify the Lie algebras to be 
su{N) (TV odd) and sdif fiT"^). Despite the weak form of convergence of 
su{N) to sdijfiT'^), solutions of the Euler equations associated to the Lie 
algebras su{N) converge, in a sense dehned further in the paper, when N 
goes to inhnity, to the solutions of the Euler equation associated with the 
Lie algebra sdiff{T‘^). Proof of this fact is the central result of this paper. 

In the fourth section we discuss numerical calculations within this ap¬ 
proach and within the Galerkin approximation. Numerical stability of this 
method in comparison to Galerkin one turns out to be signihcantly better. 
On the other hand, the numerical simulations demand signihcantly more 
calculation time than the Galerkin one. 

In the Gonclusions we discuss the hnite-mode approximation based on 
the group theoretic relations as a method of choice for solving Euler equation 
whenever numerical calculations are expected to give reliable results for long 
evolution times, and adequate computational power is available. 

We supplement the paper with an Appendix, in which we discuss ge¬ 
ometry of noncommutative tori, providing geometric interpretation of the 
approximation method. 
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2 Approximation of the algebra sdiff{T‘^) by 
su{N) algebras 

In this section we remind the formal approach to the approximation of 
sdiffiT^) and sdiff{S‘^) algebras given in [|^ and called there La-convergence. 
We remind here the resnlt on weak nniqneness (proposition 3.3 from ||13|| ). 
We also discnss behavior of generators of su{N) in the form chosen in [|I| 
when N goes to inhnity. We show that their norms diverge in this limit. 

In [|^ instead of algebras u{N), su{N) are considered their complexi- 
hcations, gl{N^(C) and sl{N^(C). It is qnite a standard to consider instead 
of algebras sdiffiT"^), sdijJ{S^), etc., their complexihcations, so the above 
procednre is fully justihed. It does not lead to any principal differences 
in comparison with our dealing with corresponding real algebras, since the 
structure constants remain the same and the relations of being isomorphic 
or non-isomorphic are preserved by complexihcation. 

Let us begin with giving necessary dehnitions of objects appearing in this 
section. First, let us remind some dehnitions and basic facts about the Lie 
algebras gl{oo), gl^{oo), and La, as they are presented in [^. 

First, gl{oo) is dehned as gl{oo) := {{aij)ij^^ : ttij G (L, all but hnite 
number of aij = 0}. This is then a matrix algebra with hnite support. As 
usual for matrix Lie algebras the commutator of matrices serves as the Lie 
bracket. 

The Lie algebra gl+{oo) is dehned in the same way with the only diherence 
that the indices of the elements of the matrices in the Lie algebra are not 
integers Z but natural numbers LV. Since Z and IN are bijective, any 
such bijection induces an isomorphism of gl{oo) and gl+{oo). However, the 
isomorphisms are not canonical, which means they depend on the choice of 
a basis in each of the Lie algebras. Since gl{oo) and gl^{oo) have only hnite 
number of elements diherent from zero, there is no problem with dehning the 
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trace of the elements in these Lie algebras. Namely, the number of non-zero 
diagonal elements is also hnite for all elements of gl{oo) and gl+{oo), and 
because of that the standard dehnition of a trace as the sum of diagonal 
elements of a matrix is unchanged for both algebras. Then the Lie algebras 
sl{oo) and sl+{oo) can be easily dehned as consisting of those elements of 
gl{oo), 5 f/+(cx)), which have trace equal to zero, and with the Lie structure 
inherited from the algebras gl{oo), gl+{oo), respectively. 

The Lie algebras sdiff{T‘^) and La, which are important in our further 
considerations, are dehned in the following way. One begins with the complex 
vector space: 

V = {Trniffl E Z X Z}, ( 1 ) 

it means a complex vector space generated by the basis T^, with indices 
rh = {mi, m2) such that mi and m2 are integers. 

The following interesting Lie algebra structures, indexed by A G IR, can 
be dehned on V: For A G iR— {0} they are dehned on the generators by: 

[Tin, Tfi]^ := sin (27rA(m x n))Tfn+n (2) 

ZnA 

and for A = 0 by: 

[Trn, Tn]° := {ffl X n)Trn+n, ( 3 ) 

where 

fh X n = min2 — m2ni for m={mi,m2), n={ni,n2). (4) 

Then, as usual, the Lie brackets are extended to the whole V by linearity. 
These Lie algebras are denoted La = {V, [•,-]'^). They are direct sums of 
subalgebras 

La = {T(o,o)} © {Lm ■ m G Z X Z — { 0 }} ( 5 ) 

where {■} denotes a vector space spanned on generators listed inside the 
hgure brackets, and the Lie structure is dehned by restriction of [•, to the 
given subspace of V. 
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The Lie algebra La is defined as the second snmmand in the above snm, 
La = {{Tfh '■ fh G Z X Z — {0}}), [•, -j^). In the case A = 0, Lo is identified 
with the complexified Lie algebra of the nonconstant divergenceless vector 
fields on T^, which is the sdijf{T'^). 

Since approximation of the Lie algebras La, sdijf, etc. by the algebras 
gl{n) or su{n) can not be realized in a strong sense of inductive or projective 
limit as it is often done for other infinite-dimensional algebras, the authors of 
[|r^ axiomatized the concept of so called La-limit of the Lie algebras. This 
notion is a generalization of intuitively defined limits from the papers [§], 

H- 

The definition consists of two axioms. The first axiom demands existence 
of a special family of real or complex Lie algebras {Qa, [•, ■]a,a £ L), where 
the index set I = IN or IR. This family should be equipped with a family 
of metrics da- Accordingly to the axiom there should exist another real 
or complex Lie algebra {Q, [•,•]) such that for every a & I there exists a 
surjective map Pa '■ G Ga and the condition 

yx,y eG if daiPa{x),Pa{y)) 0, then x = y. (6) 

should be satisfied. The second axiom concerns the definition of a Lie algebra 
structure on the La-limit. Accordingly to this axiom {G, [•, ■]), is called La- 
limit of the algebras {Ga, [•, ■]a,da,C( G I) and the family of the algebras is 
called the approximating sequence introduced by {pa, a G /) if the condition 

da{Pa[x,y],\paX,Pay]a) a=^0. ( 7 ) 

is satisfied. There (in |^) was also proved a weak uniqueness theorem, which 
states that if {Ga, [•, ‘la) da, a G /) is an approximating sequence for {G, [•, •]) 
induced by (pa, a ^ I), and G' is a linear subspace of G, which carries a Lie 
structure [■, ■]' which projects on each Ga, then {G', [•, -jO is a Lie subalgebra 
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of {Q, [•,•]) if and only if the approximating seqnence for Q is by restriction 
also approximating seqnence for Q' indnced by the restriction of the pa- 
However, this theorem abont weak nniqneness is not of very big valne, 
becanse it does not gnarantee that if we start with a given approximating 
seqnence of the Lie algebras G I) we end np with isomorphic limit 

algebras, independently of how we came to the limits. 

Let us remind how sdijJ{T‘^) is described rigorously within the La limit 
of the algebras La when A —0. One begins with the sequence of algebras 
= Li/ 7 v, i.e. the sequence of algebras La with A = L. One distinguishes 
an ideal in each algebra of the sequence: 

jN ; rn, d e ^ X (8) 


Then one dehnes the sequence of factor algebras L( 7 v) := L^/ with tta? : 

L( 7 v) the canonical projection. As the basis of each algebra in the se¬ 
quence serve 7rAr(T,n), rfi = (mi,m 2 ), 0 < mi,m 2 < N. The Lie structure 
in each algebra of the sequence is given by: 

(Ljti) , TTtv(L fj)] ^ ^ sin ^ ^ (L(m+n)mod v) • (9) 

This algebra L(Ar) is exactly identical with the algebra gl{N) with N odd. 
The approximating sequence consists then of the algebras (L(Ar), [., with 
the index set I = IN. The surjective maps Pa are identihed in this case as the 
maps ttat. The metric on each L( 7 v) is induced by dehning it on the generators 
by 

< 7rAr(T^),7rAr(Tii) >= 6 (10) 


and then extending to the whole L( 7 v) by linearity. 

It is shown in that the 2 axioms of Lc^-limit are satished for this 
sequence. Therefore, the algebra sdiff{T‘^) is an L„-hmit of the algebras 
gl{N,Cd). Nevertheless, gl{oo) which is an inductive limit of the algebras 
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gl{N,CJ), so it is the limit in a strong sense, is not isomorphic to sdiff{T‘^). 
This fact is also proved in . 


The above example shows that one should be careful with approximating 
inhnite-dimensional real algebras by hnite dimensional ones in such weak 
sense as the LQ-limit. We will demonstrate now another difficulty of this 
dehnition of a limit. The difficulty comes from the fact that the renormalized 
basis of the Lie algebra su{N) has divergent norm for N going to inhnity. 
To be more precise, let us remind that the basis of the su{N) algebra which 
leads to the mentioned -approximation of sdiff(T‘^) is dehned as: 


J. = 

where A = g = diag{l, A,..., A^“^), and 


( 11 ) 


h = 


0 1 . 

0 0 1 


0 

0 


1 . ... 0 

is a cyclic permutation matrix. Then, to obtain desired form of the structure 
constants of the algebra, a renormalization is done by multiplying all the basis 


elements by As a result, the renormalized basis elements for > oo 
become elements of inhnite norm. This is in apparent contrast to the behavior 
of the basis elements of the algebra sdiffiT"^) 


Ln = i exp {in ■ x){n 2 di — nid 2 ) 


( 12 ) 


, which are all of hnite norm. 

It seems then that the algebraic arguments on approximation of sdiff{T'^) 
by su{N) need careful treatment. 


What is also of interest, especially in the context of geometric description 
of an ideal incompressible huid, and its geometric quantization like in 
or []^ , it is a better understanding of convergence of coadjoint orbits of the 

















groups SU{N) to the coadjoint orbits of the goup SDif fiT"^). Analogous 
interesting question concerns the Casimir functions. These in the case of 
SU{N) are: trA^, trA^, ... trA^,... for an A G su{N). Their counterparts 
for the SDiffiT^) are: Jrp 2 0 j‘^n, Jrp 2 uj^^,...Jrp 2 ..., where fi is the stan¬ 

dard surface element of the torus T^. The work on this issue is in progress. 


3 Approximating Euler equation 


In this section we introduce the Euler equation for incompressible flows 
of an ideal fluid on a 2D torus. Then we remind the dehnition of the Euler 
equations connected with arbitrary Lie algebras, due to Arnold. Next, we 
specify the latter equations to the case of su{N) Lie algebras {N odd) and 
to the case of the Lie algebra sdif f{T‘^). Finally, we prove the theorem on 
convergence of solutions of the Euler equations for su{N) to the solutions of 
the Euler equation connected with the sdiff{T‘^) Lie algebra. 

As was noticed in 0 , the Euler equations can be generalized to any Lie 
algebra by 


U!. 






( 13 ) 


where are structure constants of the algebra and 5 ^ is a symmetric metric. 
For example, g is the inverse of the inertia tensor in the case of motion of a 
rigid body, for which the appropriate algebra is so(3)). The summation over 
repeating indices is understood. 

The equation for two-dimensional hydrodynamics of an ideal fluid can be 
expressed in terms of the Poisson bracket {•, •} of vorticity uj and the stream 
function as follows 

ch = a;|, (14) 

where A“^a; denotes the stream function, i.e. the solution of /S.'ip = uj. In 
the case of the flow in a rectangle with periodic boundary condition (i.e. on 
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a two-dimensional torus T^) we assume the mean vorticity equal to zero, 
fj ^2 u = 0. As noticed by Zeitlin, the symplectic diffeomorphisms of are 
generated by 


Lfi = iexp[i(niXi -f n 2 X 2 )](n 2 di — ni(? 2 ) = —iexp(in ■ x)(n x V), 

where = {n\ + n^) 7 ^ 0. These operators satisfy the Jacobi identity and 
form the Lie algebra sdif f{T‘^). For the commutators we have 


[Ln, Lm] = {n X m)Ln_^_in, where n x m = nim 2 — n 2 mi. (15) 


Equation ([T^) , when expressed in the basis (^), i.e. in terms of the Fourier 
components, {un} of cu, takes the form 


UJf. 


= E 


m X k 




~^m+k^—k J 


( 16 ) 


where = a;|, because a;(x) is real. Taking (^) as the basis for the Lie 
algebra su{N) for N = 2n + 1 one can construct a system which formally 
approximates the system (pTf) . Assuming = ^Jn-, where Jn are dehned 
in (^) we have 


TV Stt 

[Ln,Lff,] = — sin[(—)(h x m)]L 
Ztt I\ 


n+m ImodiV 


( 17 ) 


in which the summation modN in both components is understood. Thus, 
assuming the inertia tensor analogous to the one in (PIBD we may dehne the 
following Euler equations on the dual of su{N) 


sin[27r(m X fc)/A^] 

k 


^U2 


m + fe |modiV^ 


( 18 ) 


where 


N-l 


N-l 


, i = 1,2. In the literature the system (|1^) 


2 < ki,mi < 2 

is known as a hnite dimensional analog of two-dimensional hydrodynamics. 
However, no proof was known that solutions of (P!B[ ) are indeed approximating 
those of (|^). To prove that such a correspondence exists we will treat 


10 








solutions of (p!8D as infinite sequences assuming cnj;, are zero for \ ki \> 


N-l 


. Let us assume therefore that {uj'^{t)} is a solution of Eqs ([T6|) with the 
initial data {a;fc(0)}. Correspondingly, let {Coi{t)} be a solution of Eqs (|T^) 
with initial values { 0 )^( 0 )}. We will assume that the difference a;fc(0) —0)^(0) 
is small in i.e. I — ‘^fc(O) P< S. Then we will prove under some 

additional assumptions for uj that on a hnite time interval the solutions are 
also close. Let us write the equation for the difference Q(t) = ui{t) — uj{t). 
Subtracting (^) from (IT) one obtains after some manipulations 




sin(£m X k) 
efn X k 


^m+k^—k T 


sin(£m 'X k) _n t 

/ / sk"^ \S^fh+k^—k ^fh+k^^—ki 5 


( 19 ) 


where e = It is desirable to express the r.h.s. of (^) only in terms of 
variables Since 0 ;,^+^ = we have 

sin(£7Ti X k) 


= 1 

k ^ 


em X k 


—k 


^ sin(£m X A:) ^ ^ ^ ^ _n .i 

/ V ^^2 \S^fh-\-k^—k ^fh-\-k^^—k ^^in-\-k^—kj ; 

k 

Now we are ready to derive the equation for the L^-norm of f2 


1 

2 


^ ll^= I | 2 = 


2 ^ 


The last equality in the chains follows from f2(a;) being real. 
Let us observe that 

= i E f = 1 : 


21? 


dQ-f 

dt 


( 20 ) 


( 21 ) 


( 22 ) 




dt dt 

,since the sum substitution m —»• —ffi becomes 

Multiplying Eq.(^) by fl-m, performing the summation over fh one ob¬ 
tains 
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( 23 ) 


where 


^^sm{emxk) ^ ^ 

/ / / / ^^2 ™ 


m fc 


^ ^ sm{em x k) ^ ^ 

2^2-^ y2 

fh ^ 

smjemxk) 

2_^ 2-^ ^^2 ^^m+k^^—m^^—k 


m k 


ek‘^ 


'm(^) ^ ' ®m(^)^m+fc^—fc" 


and where 


i; m X k f sm{em x k) 




1 - 


em X k 


In the last two snmands repeats the term: 


( 24 ) 


( 25 ) 


4 := IIsin(em X A^)n^+jSJ_* (26) 

fh 

Changing the name of the snmmation index does not change the valne of 
Jfc. Let ns then introdnce the change from fh into —{fh + k). As a resnlt one 
gets {fh' = —{fh + k) and then we omit '). 

h = IIsin(£(-(m + fc) X 
m 

= -^sin(£(m X =-4 (27) 

fh 

Therefore Ik = 0. 

Now we will prove that the seqnence {amis)} dehned by (|2^) satishes 
S I CLm, P= S{£) 0 when £ —> 0 provided that the solntion uj{x) belongs 

to the Sobolev space[] H'"{T'^) for a > 2. We have in snch a case I P 
{l + fh?Y < oo and therefore ujfh can be represented as ujfh = where 

still Xfc I Cfc p= c < cx). Thus denoting / 3 ^{e) = a^{e){l + (rh + + 

^We know that in two-dimensional case the Euler equations of an ideal fluid have 
global solutions whose quality (differentiability properties) are related to the quality of 
initial data. 
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/c^) O'/s -^g have = Y.k Pmi^)^rh+kC-k ,which by Sch-wartz inequality gives 
us 


™ l<. EI/ 3 |(-:) P,E 


^m+k^—k 


(28) 


y /c y A: 

The last sum on the right hand side is obviously bounded because 

I Cfa+kC-k sup I Cm 1^ X! I Cfc r= Cl < oo. 

m _ 

/c k 

Therefore as it follows from 

EI««t< c. EE I 4 (e) 

m k 


(29) 


To prove that this double sum tends to zero for e —> 0, we use the integral 
criterion; we denote 

G{e) = [ dkidk 2 [ _ dmidm 2 \ P • 

J\k\>l J\in-\-k\>l 

By introducing new variables Si = {m + k) x S 2 = {rh + k) x ^, where 
k-^ = (—^ 2 , ^i) we estimate G{e) extending integration in si, S 2 on the whole 

iR2 


G(e) < [ dki dk 2 f (1 
JiR'^XD Jib? \ 


sin(£ 

1 k 


1 ^ dsi ds2 

£ 

k 1 

Si } 

' k"^ {1 + Si + slY{1 + k?Y 


where D is the unit disc on IB?, /c^ < 1. For a > 2 the integral over S 2 is 
hnite and we have 






JlR(l + sl + slr (l + sf)’-‘'“’ 

where the constant Mo- grows to inhnity for a —> 2. Therefore introducing 
z = e \k \ Si we have 


G(e) < ^ / 


dki dk2 


IB?\D fc2(l + k'^Y 


(1 _ 


^2 dz 
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2 

Noting that | 1 — |< min (^,2) we can estimate the integral over dz 

splitting it in two parts as follows 


r°°, sin z, 9 dz 

(1- — 


< 


/o 


(1 + ~ 36 


(■ fc£|“ 

poo 

/ z^dz + 

4 /_ \ek z^-^^dz 

Jo 

J\k£\°^ 


1 


1 ke r" 4 

-r ek 


1 2 - cr 1 ' ' 


\2<j-l+2a{2-a) 


Taking a = we make both exponents equal to 7^ and hnally one 


2(7-1 


2cr+3 


2(t+3 


obtains 






k 


\la-i 


1126 cr — 2 i JiB?\D (1 + k?) 

.1 U2 


dk^ dk'^, 


where r = 8^^. For a > 2, the integral over k^, k^ is hnite, and thus 
G{e) ^ 0 for £ —>• 0. 

In this way we have proved that for a > 2 the hrst term on the r.h.s. of 
Eq(P^ is vanishing when N ^ oo (e = 2 tt/N). It converges to zero as 
i.e. as where r depends on a. Indeed, we have 


1/2 


'^afn{e)n_fh |<| ^ I (511 n =1^1 


Now let us estimate the second term on the r.h.s of Denoting ft = rfi + k 
we have 

sin(e/i X k) p, X k Cp 

j ^2 ^fh+k^^—k^^—fh I 2 , I 

fh k ^ 


^-k 1 ^ 51 

fi,k 


e{pxk) F (i + ;,2)(7/2 


511 c/i (1 + /i^) 

p.,k 


n_-,n-,_p\<K^ II D 


where K^j = J2p I (1 + |. From Schwartz inequality is hnite for 

cr > 2. 

Now let us come back to Eq.(^^. Dividing both sides by || fl{t) || and 
integrating with respect to time we come to 

pt /jf 

II 11 = 11 ^(0) II + / II Q/,N II 5^«m(£)hlm(^) + 

Jo II D(t) II rh 

p dt sin(em x k) n n 

Jo 

" V / II ra^k 
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Applying the estimations which we had derived one arrives at 

II ^{t) ||< {|| f^(0) II +tG^/\e)} + KJ Q{t) \\dt, 

which by Gronwall inequality implies that on the time interval (0, T) we have 

II ^{t) ||< {|| G(0) II 0 < t < T. 

For iV —> cx) the difference ^2(0) between a;(0) and d;(0) tends to zero, at 
least as ~ || f2(0) ||~ Q, similarly, G(e) ~ N~^ 0, 

and therefore also || Q{t) ||—*> 0, which proves that the solution ui of Euler 
equation on su{N) approaches the corresponding hydrodynamic solution oj. 

In a similar way one can prove a little more general theorem which sum¬ 
marizes our considerations. 

Theorem If a;(t) G C'^([0,T], is a solution of the Euler equation of 

hydrodynamics on the torus then corresponding solutions (i.e. such 
that d;,^(0) = cc;,^(0) for < rrii < of the su{N) Euler equations 

are approaching oj in the norm of for every k satisfying 0 <k < a —2. 

4 Numerical simulations 

The theorem proved in the last section shows that for large N the difference 
between solutions of the Euler equation on su{N) and the Galerkin approxi¬ 
mation of the Euler equation for ideal incompressible fluid is negligibly small. 
Numerical schemes usually applied to solve Euler equation are based on the 
Galerkin approximation. Numerical calculations were performed to test if the 
group-theoretic approximation could be applied to create a better numerical 
scheme for solving the Euler equation. There seems to be one important ad¬ 
vantage of the group-theoretic hnite-mode approach over the Galerkin one. 
Namely, the former possesses a big number of conserved quantities, which 
grows with N. The Galerkin approximation does not have this property. As a 
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result of this difference one can expect higher stability of the group-theoretic 
approximation for long times. 

In the computer simulations we adopted similar initial conditions as 
Dowker and Wolski |^. This means we took uj(x,t = 0) ~ N‘^Sx,xo ~ 1- 
Strictly speaking, such initial conditions are not satisfying the assumptions 
of our theorem. However, in both approximations, it means the Galerkin and 
the group- theoretic ones, these initial conditions lead to numerical solutions 
quite well approximated by the discrete ones. Dowker and Wolski ^ were 
performing their simulations at the level of stream functions, we were doing 
that at the level of vorticities. 

Calculations in Q] were done on a Hewlett-Packard workstation. We 
agree with that it is very difficult to go to Ns big enough for a convinc¬ 
ing numerical illustration of the convergence proved in the theorem. The 
complexity of the system of ODEs strongly depends on N. In our case the 
number of equations grows linearly with N, the number of quadratic terms 
in every equation grows with The N is limited not only by the execution 
time, but also, and in fact mainly, by the memory requirements. This is why 
we hope that our method will be much more effective in the future — we 
witness every year signihcant progress in memory availability. 

It was not our aim at this point to make calculations for the biggest 
possible N, but rather to observe the numerical stability of both methods. 
This is why calculations were done on a PC (Pentium H, 700MHz, ICB 
RAM), in Mathematica 4.1. We restricted ourselves to small Ns (relatively 
easy simulations). The effect is sufficiently strong to be observed at = 
11,15,25. It is expected from theoretical considerations to be much more 
signihcant for big N. 

We observed individual Fourier modes as functions of time for both ap¬ 
proximations. We were testing time after which the symmetry a;| = 
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(for our choice of initial conditions (real uij.) it was uij. = was breaking. 
One can easily observe big difference in times of destabilization between the 
two approximations. For iV = 11 and iV = 15 numerical calculations based 
on the Galerkin method destabilize after the amount of time approximately 
equal to 4.5 computer units, while the group-theoretic method presented 
stable behavior till 20 computer time units. For N = 25 the difference was 
signihcantly bigger. Namely, the group-theoretic method destabilized after 
about 18 units of time, while the Galerkin method destabilized after 1.5 
units of time. In this case the group-theoretic method was more stable than 
the Galerkin one by a factor of 12. This difference was signihcant, and we 
are convinced that with bigger N and longer simulation times the difference 
would increase. We used in both cases exactly the same initial conditions, 
numerical methods, numerical accuracy and all other parameters. 

As an example, the results of the numerical simulations for the represen¬ 
tative cases of N = 11,15, 25 for a pair of modes each in the Galerkin and 
in the group-theoretic approximations are shown in Figs. 1-3. In the graphs 
is presented the time dependence of the amplitudes of three particular pairs 
of modes with the numbers k and —k. They are representative for all the 
pairs of modes. We were looking for the instant of time when the two modes 
begin to behave differently showing that reliability of the simulations is lost. 

5 Conclusions 

In this paper we studied the approximate ideal fluid description based on 
the observation concerning the convergence of the structure constants of the 
su{N) {N odd) algebras to the structnre constants of the algebra sdiffiT"^), 
and on similarity of the Euler equations dehned on the dual algebras of su{N) 
and sdijfiT^), respectively. The hypothesis stated in the literatnre that the 
Euler equations for su{N) converge to those for sdiffiT"^) was carefully ex- 
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amined. 


The results are: 

In addition to the fact known from the literature that the algebras 
su{oo) and sdiffiT"^) are not isomorphic, we pointed another difficulty of 
the approximation. Namely, we have shown that in the limit N ^ oo renor¬ 
malized elements of the basis of su{N), used to demonstrate similarity of 
the strucure constants of lim su{N) and of sdif f{T‘^), become elements of 

N^oo 


infinite norm. 

We proved the statement about the convergence of the solutions of the Eu¬ 
ler equations for su{N) to the solutions of the Euler equations for sdif fiT"^). 
This way the functional analytic part of the approximation eventually found 
a complete proof. 

We performed extensive numerical simulations based on the group-theoretic 
and the Galerkin finite-mode approximations. They showed that the group- 
theoretic method is significantly more reliable for long evolution times than 
the Galerkin one. When the same numerical conditions are applied, the 
Galerkin method leads to instabilities after evolution time at least a few 
times shorter than the group-theoretic one. This difference grows with the 
number of modes used in symulation. For N=25 group-theoretic numerical 
solutions were 12 times more stable than the Galerkin ones. We expect the 
difference to be sharper for bigger numbers of modes involved. On the other 
hand, the Galerkin method is significantly less machine time and memory 
consuming than the group-theoretic one. Galerkin method is appropriate 
for short evolution time calculations, performed with limited computational 
power. With given machine power at hands one can obtain solutions for short 
evolution time using the Galerkin approximation for bigger N than using the 
group-theoretic one. 

The long time scale solutions of the Euler equations are used in meteo- 
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rology, in weather prediction, as well as in studying oceanic flows, to name 
just a few applications. The group-theoretic method would allow to improve 
long time scale predictions reliability in the mentioned areas. It currently re¬ 
quires powerful workstations for implementation. With dramatic growth of 
computational power and the memory size in recent years, in a few years the 
group-theoretic method should be easier to apply than it is at the moment. 

6 Appendix 

Noncommutative tori are one of the simplest examples of noncommutative 
spaces, which are a subject of study in noncommutative geometry. In this 
appendix we present material on noncommutative tori, which is relevant for 
the finite-mode approximation discussed in the paper. 

The paradigm of noncommutative geometry grew up from the Gelfand- 
Naimark theorem, accordingly to which the category of Hausdorff, locally 
compact topological spaces is dual to the category of commutative G*-algebras. 
The commutative G^-algebras are algebras of continuous complex-valued 
functions (with compact support in the case of noncompact spaces) on the 
spaces. The G*-algebras of continuous functions on compact spaces are ad¬ 
ditionally unital, with the unit given by the constant function equal to one. 
Constant functions cannot have compact support on noncompact spaces, 
therefore the corresponding C'*-algebras are nonunital. Thanks to the Gelfand- 
Naimark theorem, the entire topological information on the spaces is coded 
in the algebraic properties of commutative C'*-algebras, to which they are 
dual. 

Relaxing the commutativity property on the G^-algebras locates us then 
in the held of noncommutative topology. The subject of its study is not some 
kind of noncommutative spaces, but rather algebras of (continuous) functions 
on them. These ’’functions” are actually operators, and the G^-algebras 
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are algebras of operators acting in appropriate Hilbert spaces, in general 
infinite dimensional. This statement is the content of a second theorem 
by Gelfand and Naimark. (On Gelfand-Naimark theorems one can consult 
for example the recent book |^). The fact that coordinates of the virtual 
noncommutative spaces are not commuting, and the possibility to realize the 
algebras of ’’functions” on them as algebras of operators in a Hilbert space, 
inspired the creators and developers of this area of mathematics to adopt 
physicists’ terminology and to speak about quantum geometry or quantum 
topology (the latter name in the case of strictly limiting to the G^-algebras 
seems more appropriate) and about quantum spaces. 

In contrast to topologists, geometers usually prefer to speak about smooth 
manifolds rather than topological spaces. Therefore, instead of G^-algebras 
noncommutative geometry studies ^-algebras that are dense in appropriate 
G^-algebras exactly in the manner in which algebras of smooth complex val¬ 
ued functions on manifolds are dense in corresponding algebras of continuous 
functions. The conditions lead to the notion of a pre-C'*-algebra. 

The noncommutative torus is a quantum space, which might be described 
using either C'*-algebraic or ^-algebraic language, with appropriate limits in 
the latter case locating us in the C'*-algebraic framework. Let us shortly 
review, following [^, the ^-algebraic approach to the noncommutative 
torus. The algebra Ae of smooth ’’functions” on the noncommutative torus 
Tq is the unital ^-algebra generated by 2 indeterminates Ui, U 2 with the 
relation 

f/if /2 = (30) 

The * structure is defined by imposing unitarity on Hi, 1 / 2 , by Ul = 

U 2 = a generic element a G ^0 is then of the form 


a = 




n 


( 31 ) 
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, where a^n is a complex valued Schwarz function on i.e. a sequence 


of complex numbers {amn € (C : {m,n) G which decreases rapidly at 

infinity (\/k G IV ||a||fc := swiJ^lamnKl + + |^|)^ < oo ). 

In this paper we are interested in the case, in which 6 is rational, 6 = 
M/N, where M and N are positive integers taken to be relatively prime. 
Therefore we present further only the rational 9 case. In this case Am/n 
is the algebra of smooth sections of a twisted matrix bundle over T^. Let 
us briefly describe the construction leading to the twisted bundle following 
The bundle is denoted Aiq — >• T^, with q = . The fibers of the 

bundle are isomorphic with the algebra Mat{N,^J) of all iV x iV complex¬ 
valued matrices. Its transition functions have values in GL[N,^). They act 
on Mat{N,G) by conjugation. 

First, we consider a trivial bundle (we identify here and below with 

X S^) 

X X Mat{N,(D) x (32) 


where pri 2 is the natural projection on the first two members. The algebra 
of smooth sections of this bundle is 


X S\Mat{N,([!)) = C°°{S^ x S\(D) ® Mat{N,(D) (33) 

The algebra Mat{N,G) may be understood as the unique algebra gener¬ 
ated by 2 unitary elements Uq and Vq such that 

UoVo = qVoUo, = 1 (34) 

The matrices Uq and Vq may be represented in the form, which is very 
closely related to the matrices g and h introduced in the section 2. Here 
however, the parameter q differs from the parameter A, since q = but 

A = . The algebra of smooth sections of the bundle (^) is generated 

by unitary commuting elements u, u, and by unitaries Uq, Vq generating 
GL{N,(C). 
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In order to define the twisting of this bundle, the following actions should 
be considered: the action of X on x {u, v) i—>• {q'^u, q'^v), 

m,n & Zn, and the action of on Mat{N,^), A i—> UqVq^AV^Uq"', 
m,n E 

Then, we impose the condition that the following diagram is commutative 

^2 

X X Mat{N,(U) -A Mg 

Wl,2 i [Pq (35) 

772 

S^xS^ ^ S^xS^ 

There is a vector bundle over x , £q ^ x such that Mq = 
End{£q). This fact is closely related to the Morita equivalence of the non- 
commutative torus with the commutative one. Namely, Morita equivalence 
of two C^-algebras or two pre-C^-algebras means that the algebras have a 
common module, on which one of them acts from the left and the other 
from the right. This is exactly the case for a noncommutative torus with 
rational ’’deformation parameter” and the standard commutative torus. We 
may treat the space of (smooth) sections of the bundle £q ^ x as the 
left module for the algebra Mq, since the latter consists of endomorphisms 
of the former. Simultaneously, since £q ^ x is a bundle over the 
standard torus, the pre-C^-algebra of smooth functions on the torus acts on 
the space of smooth sections of £q from the right. This is the way in which 
the rational noncommutative tori are related via Morita equivalence with the 
commutative torus. 

However, the Morita equivalence of the two kind of spaces, although very 
important, is not directly responsible for the convergence effects we studied 
in the body of this paper. In the geometric picture the limiting procedure 
may be described in the following way: The algebra Mat{N,^), which is 
isomorphic with the hbers of the bundle Mq ^ x may be undestood 
as the Lie algebra gl{N,^). From this moment on one should just remind 
the construction described in the body of the paper. 
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Figure captions 


Fig 1.1 Group - theoretic method, N=ll, the modes compared: 

i) mode (—5, 0) and mode (5, 0) 

ii) mode (—3, 2) and mode (3, —2) 

iii) mode (—3, 0) and mode (3, 0) 

The simulation was still stable at 20 units of time. After this time the 
k and — k modes started to differ from each other, which means that 
the solution is no longer stable at this point. 

Fig 1.2 Galerkin method, N=ll, the modes compared: 

i mode (—5, 0) and mode (5, 0) 

ii mode (—3, 2) and mode (3, —2) 

iii mode (—3, 0) and mode (3, 0) 

The simulation was stable till 4.8 units of time. After this time the k 
and —k modes started to differ from each other, which means that the 
solution is no longer stable at this point. 

Fig 2.1 Group - theoretic method, N=15, the modes compared: 

i) mode (—7, 0) and mode (7, 0) 

ii) mode (—1,2) and mode (1, —2) 

iii) mode (—7, 5) and mode (7, —5) 

The simulation was stable till 23 units of time. After this time the k 
and —k modes started to differ from each other, which means that the 
solution is no longer stable at this point. 

Fig 2.2 Galerkin method, N=15, the modes compared: 
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i) mode (—7, 0) and mode (7, 0) 

ii) mode (—1,2) and mode (1, —2) 

iii) mode (—7, 5) and mode (7, —5) 

The simulation was stable till 2.7 units of time. 

Fig 3.1 Group-theoretic method, N=25, the modes compared: 

i) mode (—12,2) and mode (12, —2) 

ii) mode (—6, 2) and mode (6, —2) 

iii) mode (—5, 2) and mode (5, —2) The simulation was stable till 18 
units of time. After this time the k and — k modes started to 
differ from each other, which means that the solution is no longer 
stable at this point. 

Fig 3.2 Galerkin method, N=25, the modes compared: 

i) mode (—12,2) and mode (12, —2) 

ii) mode (—6, 2) and mode (6, —2) 

iii) mode (—5, 2) and mode (5, —2) 

The simulation was stable till 1.5 units of time. After this time the k 
and —k modes started to differ from each other, which means that the 
solution is no longer stable at this point. 


27 














o Ms ('JfO)jCl,o} 












v^oXeS n,l) 










^ j «^octes (- 1-jO) , C'?',0') 








7.5 





vnodes 









5 











10 




(\j' -S. /| VW> b y ) , (s>a^t''^^'‘ 






wao 






6-6 [2.) iCCi-2.') J cSaleAr-kt^K 






fs[-2.5j C ^1^); 




-8 


svioctes (^-'12,2)^(/12^-2.') ^ 








-i 





